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Abstract 

Optimizing the energy efficiency (EE) for the MIMO broadcasting channels (BC) is addressed in 
this paper, taking into account the transmit independent power which is related to the active transmit 
antenna number A new optimization framework is proposed, in which transmit covariance optimization 
under fixed active transmit antenna sets is first performed and active transmit antenna selection (ATAS) 
is utihzed then. To optimize the EE under a fixed transmit antenna set, we propose an energy efficient 
iterative waterfilling scheme according to the block-coordinate ascent algorithm, through transforming 
the problem into a concave fractional optimization via uplink-downlink duality. It is proved that the 
proposed scheme converges to the global optimality. After that, ATAS is employed to determine the 
active transmit antenna set and to turn off the rest inactive antennas. ATAS can balance the active 
transmit antenna number related EE gain with higher capacity gain and the EE loss with more transmit 
independent power wasting. During the ATAS, the optimal exhaust search and norm based successive 
selection schemes are borrowed. Through simulation results, we discuss the effect of different parameters 
on the EE of the MIMO BC. 
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I. Introduction 

Wireless communication turns to the era of green. This is not only because of the exponential 
traffic growth with the popularity of the smart phone but also the limited energy source with 
ever higher prices. Energy efficiency (EE), as a result, becomes one of the major topics in the 
research of wireless communications [2J and plenty of research projects either government funded 
or industrial funded start to investigate the energy efficient solutions for the wireless network 
as well as the sustainable future for the wireless communications. Meanwhile, multiple input 
multiple output (MIMO), especially downlink multiuser MIMO (also called MIMO broadcasting 
channels, BC), is becoming the key technology in the next generation cellular networks due to 
its significant spectral efficiency (SE) performance. Although the SE for the MIMO BC is well 
known, the EE still remains an open issue. Thus, studying the EE of the MIMO BC is important 
for the design of the future green wireless networks. 

The capacity of the MIMO BC has been widely studied in the literatures, e.g. [I3]|-[l5]| and 
it is well known that dirty paper coding (DPC) can achieve the capacity region of the MIMO 
BC. Through transforming the capacity of MIMO BC into the dual MIMO multi-access channel 
(MAC) dH, the capacity can be efficiently optimized through the convex optimization and a 
well structural solution is called iterative waterfilling [|51, O. However, the EE of the MIMO 
BC has been rarely studied. The EE is in general defined as the capacity divided by the power 
consumption, which denotes the delivered bits per-unit energy measured in bits per- Joule. For the 
MIMO BC, when the optimization objective transfers from the capacity to the EE, there is a two- 
fold impact. On the one hand, the optimization varies from a convex optimization problem to a 
fractional programming problem [7J, when the uplink-downlink duality is applied. The previous 
spectral efficient schemes are not practical any longer. On the other hand, activating all available 
transmit antennas and radio frequency (RE) chains for transmission is not always energy efficient 
although spectral efficient, as more transmit antennas correspond to higher transmit-independent 
power, which induces a tradeoff between the EE gain with higher multiplexing/diversity and the 
EE loss with more transmit independent power consumption wasting. 

The EE of the MIMO BC is studied in this paper, where a practical power model with transmit 
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independent power |[8l, BUl is taken into account. The transmit independent power contains the 
cost of signal processing, circuit power, etc. at the BS, and is highly related to the active 
transmit antenna number. A new optimization framework with transmit covariance optimization 
and active transmit antenna selection (ATAS) is proposed to maximize the EE of the MIMO BC. 
At first, we find that the EE optimization problem under fixed active transmit antenna set is a 
concave fractional programming, and propose an efficient energy efficient iterative waterfilling 
scheme to obtain the optimal covariances, which is proven to be globally optimal. After that, 
exhaust search and norm-based successive selection are borrowed in ATAS to determine the 
active transmit antenna set and turn off the inactive ones. 

A. Contributions 

We observe that the EE is affected by both the transmit covariances (total transmit power) and 
the active transmit antenna sets. On one hand, the active transmit antenna sets are related to the 
transmit antenna selection diversity and thus affect the capacity. On the other hand, the size of 
the active transmit antenna sets is equal to the active transmit antenna number, which is related 
to the transmit independent power. Based on this observation, a new optimization paradigm with 
transmit covariances optimization and ATAS is proposed. 

Under fixed active transmit antenna sets, we derive the optimal energy efficient transmit 
covariances at first. Employing the famous uplink-downlink duality, the nonconcave MIMO BC 
capacity is transformed into the dual concave MIMO MAC. Correspondingly, the EE becomes a 
quasi-concave function. After that, we propose a novel well structured energy efficient iterative 
waterfilling scheme based on the block-coordinate ascent algorithm to solve the EE optimization 
problem efficiently. During each iteration, the transmit covariance matrices optimization is formu- 
lated as a concave fractional program, which is solved through relating it to a parametric concave 
program and applying the Karush-Kuhn-Tucker (KKT) optimality conditions. Interestingly, the 
solution of each iteration has a feature of waterfilling. We prove the convergence of the proposed 
scheme and validate it through simulations. 

A novel ATAS procedure is proposed to compromise between active transmit antenna number 
related EE gain with increasing multiplexing/diversity and EE loss with increasing transmit 
independent power consumption. Exhaust search and norm based successive selection schemes 
are borrowed to choose the energy efficient active transmit antenna sets. A unique feature of the 
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ATAS is that the inactive antennas should be switched off to save power, e.g. employing micro- 
sleep ifTOl or discontinuous transmission (DTX) [fTT|. Furthermore, the invisible channel state 
information at transmitter (CSIT) problem during the implementation is discussed. Simulation 
results give us insights about the effect of different system parameters on the EE. 

B. Related Works 

There are a lot of literatures discussing the EE of the point to point MIMO channels with 
transmit covariances optimization without antenna selection [fT2 | - [fT7l . The point to point MIMO 
channels can always be separated into parallel sub-channels through singular value decomposition 
(SVD) or after detection. In this case, only power allocation across the sub-channels needs to be 
optimized to maximize the EE lfT3l - lfT5l . As the sub-channels are parallel, the solution is similar 
with the energy efficient power allocation in OFDM systems ifTSl . |fT9l . The optimization for 
point to point MIMO channels is not applicable for the MIMO BC, as the MIMO BC cannot 
be simply transformed into parallel sub-channelil]. There are few literatures discussing the EE 
for the MIMO BC. To the best of the authors' knowledge, only [201 and our previous work [iBl 
addressed the EE of the MIMO BC, but they both assumed linear precoding design and equal 
transmit power allocation for simplification. These assumptions make the both works far away 
from the optimal solution. To optimize the EE of the MIMO BC, deriving the well structured 
transmit covariances is a challenge. 

Antenna selection is a widely discussed technology in spectral efficient MIMO systems, both in 
transmitter side and the receiver side, e.g. in [[2T | - ll25l . However, the spectral efficient transmitter 
antenna selection [|2T1l - [|23l is always performed to choose the active antennas when the number 
of radio frequency (RF) chains is larger than the number of the antennas. Meanwhile, the receive 
antenna selection [[24|. [|25l is always performed for the linear precoding schemes to approach 
the asymptotic optimal performance. In our scenario, as DPC is employed, the receive antenna 
selection is unnecessary. Moreover, as we consider the case when the transmit antenna number 
is equal to the RF chain, the purpose of the transmit antenna selection is saving power through 
turning off the inactive RF chains. Although the selection procedure with exhaust search and 

'Although after zero-forcing (ZF) precoding, for example, the MIMO BC can be separated into parallel sub-channels, the ZF 
scheme is far away from the optimal solution l3l . 
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Norm based successive scheme is borrowed from the previous literatures, it is in some sense 
different from the spectral efficient ones. Especially, there exists another challenge of invisible 
CSIT problem, which is also discussed in this paper. 

C. Organization and Notation 

The rest of this paper is organized as follows. Section |n] introduces the system model and the 
problem formulation. Section UlI] proposes the transmit covariances optimization to maximize the 
EE under fixed active transmit antenna set. Section |IV] proposes the energy efficient ATAS and 
discusses the implementation issue in the realistic systems. The simulation results and discussions 
are given in Section |Vl Finally, section |VI] concludes this paper. 

Regarding the notation, bold face letters refer to vectors (lower case) or matrices (upper case). 
The superscript H and T represent the conjugate transpose and transpose operation, respectively. 
Tr(-) denotes the trace of the matrix. 

IL System Model and Problem Formulation 

The system consists of a single BS with M antennas and K users each with N antennas, 
which is shown in Fig. [T] We assume that the number of RF chains is equal to the number 
of antenna]^. Denote the channel matrix from the BS to all users as H G C^^^*^ with H = 
[Hf , Hi, . . . , H^]^, where Hj G C^^^"^ is the channel matrix from the BS to the ith user. As 
active antenna number at the BS has significant affect on the EE, selecting the active transmit 
antennas is important. Consider that the selected transmit antenna set is T C {1, . . . , M} with 
active transmit antenna number Ma = |T|, and denote the channel matrix from the BS's active 
transmit antennas to the users as H7- G C^-^^*^", with H7- = [H^^,H^25 • • • ^'^tkV^ where 
iiT,i £ C^^^" is the channel matrix from the BS's active transmit antennas to the ith user. 

The downlink channel can be denoted as 

Yi = Hr,iX + nj,i = (1) 

^The results here can be extended to the general case with different antenna number at each user and are also applicable to 
the multi-cell scenario with BS cooperation. Moreover, if the number of RF chains is smaller than the antennas, our results can 
be simply extended after some modifications of the ATAS. 



Februaiy 17, 2012 



DRAFT 



6 



where x G C is the transmitted signal on the downlink, Xj G C is the transmitted signal 
of user i on the uplink, n,; G C^^^ and n G C^^"^^ are the independent Gaussian noise with 
each entry CA^(0,cr^). Frequency flat fading channels with bandwidth W is considered and the 
channel state information (CSI) is assumed to be perfectly known at the transmitter and receivers. 

About the power model, as BSs take the main power consumption in the cellular networks, 
the users' consumed power is not considered here. The power radiated to the environment for 
signal transmission is only a portion of its total power consumption [9J, so the practical transmit 
independent power including circuit power, signal processing power, cooling loss etc. at the BS 
should be taken into account. The transmit independent power for the BS deployed with multiple 
antennas is mainly related to the active transmit antenna number. Thus, as a good approximation 
of the practical power model, we consider a general model given as 

Ptot.l = f{P, Ma), (2) 

which is mainly related to the transmit power and antenna number. We can assume that /(P, Ma) 
is monotonously increasing as a function of P and Ma, respectively and based on li26l , we assume 
that f{P,Ms) is affine or convex as a function of P. Motivated by |[8l, lIH, more specifically, 
we consider an affine power consumption model, which can be denoted as 

P 

-Ptotal = 1- MaPdyn + -Psta, (3) 

V 

where r] denotes the power amplifier (PA) efficiency; MaPdyn denotes the dynamic power 
consumption proportional to the number of radio frequency (RF) chains, e.g. circuit power of RF 
chains which is always proportional to M^; and P^ta accounts for the static power independent 
of both Ma and P which includes power consumption of the baseband processing, battery unit 
etc.. MaPdyn + Psta is referred to as the transmit independent power. 

Note that although the optimization procedure is performed based on the affine model, the 
idea can be simply extended to other convex power model case, e.g. considering the non-ideal 
PA efficiency like [|27ll or rate dependent P^ta like [l28ll . Meanwhile, note that we omit the effect 
of the complexity of the algorithms related signal processing power in the power model, as it is 
practical that the Psta and Pdyn would be significantly larger than the signal processing power 
of the algorithms. 
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A. Problem Formulation 

The EE is defined as the sum capacity divided by the total power consumption. Considering 

the maximum transmit power constraint at the BS and the minimum sum capacity constraint, 

we can define the optimization problem as 

Cbc (Hr,i, • • • , Hr,K, -P) 
max F = 5 

r.PMa I + M.Fdyn + ^ta 

S.t.O < F < P„,ax ^^-^ 

Cbc (Hr,i, • • • , Hr.i^, P) > Cmm- 

In the above problem, the optimization should be performed via optimizing P, Ma and T. 

Look at the three parameters. P affects the input transmit covariances directly. As Ma = |T|, 

the chosen of T would affect not only the channel matrices, but the dynamic power consumption. 

Therefore, the problem should be solved by optimizing the transmit covariances and selecting 

the optimal transmit antenna set. Rewrite the optimization problem as 

^ C'bc (Hr,i, • • • , H-T-i^, P) 
maxf = 5 ; — ; 

r,P f + |r|Pdyn + Psta 

S.t.O <P< P^ax 

Cbc (Hr,i, • • • , iiT,K: P) > Cmin- 
Interestingly, the structure of the above optimization problem is similar to the capacity max- 
imization of the MIMO BC with linear precoders, e.g. Il24l . Il25l . ||29ll , where power allocation 
with precoding and user/receive antenna selection are performed as two steps to find the optimal 
solution. Motivated by this idea, to get the optimal solution, the following two-step structure 
with transmit covariances optimization and the ATAS should be employed. At first, under fixed 
set T, optimal transmit covariances can be derived. And then we should choose the most energy 
efficient T. 

It is worthwhile to note that the user/receive antenna selection is unnecessary here. As during 
the optimization of transmit covariances under DPC, the selected users/receive antennas are 
correspondingly determined. If we consider the linear precoding schemes, e.g. zero-forcing 
precoding, joint transmit/receive antenna selection and precoding with power allocation should 
be taken into account, an example can be found in our previous work BOil . 
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B. Achievable Sum Capacity and Energy Efficiency 

The sum capacity of the MIMO BC is achieved by DPC, which can be denoted as follows 
jSl given a total transmit power P. 



I + ^Hr,iSr,iHr,i 



+ 



Cbc (Hr,i, . . . , iir,K, P) = max W log 

{Sr,.}|l^:Sr,.>0,E£i Tr(Sr,0<^' 

|l+ ^Hr,2 (5]r,i + ^r,2) Hf 2I |l + ^Hr,i^ (Sr,i H hSr,x)H 

M^log-! if 7— -— — , I H hlVlog 



I + ^Hr,2 (Sr,i) 2I |l + ^Hr,A' {^r,i + ■■■ + ^t,k-i) ii'r,K\ 



(6) 

where the optimization is performed to choose the optimal downlink transmit covariance matrices 



Taking ^ into the optimization problem ([5]). As ([61) is nonconcave, even optimizing ^ under 
fixed T is nontrivial. Fortunately, motivated by ll26l . [|3T1| . we find out that the following 
property. If the numerator (sum capacity) can be transformed into a convex function, the EE can 
be formulated as a quasiconcave function, because the denominator (total power consumption) is 
affine (also holds for the convex case). Based on this observation and the famous uplink-downlink 
duality, we can transform the EE into a quasiconcave function under fixed T. 

The dual uplink channel is denoted as 

K 

Ymac = E _.Xj + n. (7) 

i=l 

Applying the uplink-downlink duality [4J, the MIMO BC sum capacity Q is equal to the concave 
sum capacity of the MIMO MAC with sum transmit power constraint, which can be denoted as 



1 ^ 



CuAc {Ht,!, • • • , iir,K, P) = max W log 

{Qr,«}'l^:Qr,«>0,E£iTr(Qr,>)<P " i=i 

(8) 

where the uplink transmit covariance matrices Qr,i G C^^^, z = 1, . . . , K need to be optimized. 
Thus, the dual MAC optimal EE should be rewritten as 

c /XT XT ^ • • • ' Hf,A'7 P) 
^MAC lHr,i, • • • , -Hr,i^) = max p — , 

/ ^ + I ' K dyn + -Tcon 

(9) 

S.t.O < P < Pmax ^ ' 

CmAC (Hr,l, • • • , ^T,K, P) > C'min- 

According to the duality and the mapping between Sf i and Qt- ,; [4J, optimal Sr,*- * = 1, • • • , -f^ 
and P can be obtained if we can get the optimal Qr,i, i = I, ■ ■ ■ , K and P in (|9l). Furthermore, 
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^ Tr (Qj) = P is always required for optimizing ([8]) Q, dH) can be simplified and rewritten as 



lyiog 

^MAc (Hr,i; • • • > Hr,E-) = max rnax 



K 

i=l 



K ' (10) 

S.t.O < ETr(Qr,i) < ^max 
i=l 

CmAC ^Ht-^i, . . . , 'H.T,K, ^ Tr (Qt-.i)^ > C'min- 

Finally, the optimization of (|5]) is transformed into optimizing (flOl ). If we can obtain the 
optimal Qr,i5^ = 1, . . . , A' for (flOl) . the optimal P can be decided correspondingly based on 

K 

Tr (Qr,j) = P- As (flOl) is a concave fractional programming problem under fixed T, we 

i=l 

will separate the optimization into two steps. First, we can find the optimal transmit covariances 
under fixed T and then applying the ATAS to choose the optimal T. 



III. EE Optimization Under Fixed Transmit Antenna Set 

In this section, we will solve the optimal energy efficient transmit covariances under fixed T. 

Let us look at (flOl) with fixed T again. Since the numerator is concave and the denominator is 
affine, (flOl) is a quasiconcave optimization, which can be solved through the bisection method or 
interior-point methods ll32l . However, the numerical methods would be still too complex when 
the user number becomes significantly large. Motivated by 0, 0, an energy efficient iterative 
waterfiUing is proposed in the next section to solve it more efficiently. For ease of description, 
we omit T in the subscript in this section. 

We will first consider the solution of (fTOl) without the minimum capacity and maximum 
capacity constraints in subsection IIII-Al and then we will extend it to the general constrained 
case in subsection IIII-BI 

A. Unconstrained Case 

1) Motivation: As the EE is distinct from the capacity, the spectral efficient iterative water- 
filling |l5l is not applicable for the EE any longer. Nevertheless, we notice that the basic idea of 
the spectral efficient iterative algorithms are based on the block-coordinate ascent algorithm [[33l 
Sec. 2.7]. That is to say, if we can write the EE as the similar structure with the block-coordinate 
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ascent algorithm and then prove it satisfies the condition of |l33l Sec. 2.7], we can also obtain 
an iterative solution of the problem (flOl) . 

For ease of description, we define the following function g (■) at first. 



fif (Qi, . . . , Qa-) 







W\og 


z — 1 



(11) 



+MPdyn + ftta 



For the block-coordinate ascent algorithm, given the current iterate Q^'^^ = ^Q^^'', . . . , Q^'')' 
the next iterate Q('^+^) = ^Q^'^^^'', . . . , Q^^^^j can be generated as 



Q,- = arg max g Q 
' Q^■.Q^>o \^ 



) • • • ) ) ^J) ^j + l5 ■ ■ ■ 5 



(12) 



However, to apply the iterative algorithm efficiently, there are conditions need to be satisfied. For 
one thing, the solution of (fT2l) should be uniquely attained tSS] Proposition 2.7.1]. For another, 
the solution should be simple and easy to employ. 

Very fortunately, the two conditions both fulfill and the solution can be obtained following an 
energy efficient waterfiUing feature. We are interested to show it in the next subsection. 

2) Energy Efficient WaterfiUing: Based on [|5l, [|6]|, it is fulfilled that 

K 

E 

i=l 



log 



log 



+ log 



I + ^ E Hf Q,H, 

-1/2 



-1/2 



I + cr'I + E Hf Q,H, 



xHfQ,H,|cr2l+^HfQ,H, 



log|Z,|+log|l + GfQ,G,|, 



^1/2 



(13) 



where = I + ^ ^ Hf Q^H^ and = a^I + ^ Hf Q^H^- . By denoting 



ETr(Q,) 



6, = W^log|Z, 
and substituting (fT3l) into ([TT]) we have that 



fir(Ql,...,Qx) 



6, + Wlog|l + GfQ,GJ 



Tr(Q,) 

7] 



+ ai 



(14) 
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Therefore, we can redefine the problem (fT2l) by removing the iteration number as to 

^ bi + W\og\l + G^QiGi\ 
maximize 5f(Qi,...,Qi_i,Qi,Qi+i,...,Q/<) = —r^ (15) 

by treating Qi, . . . , Qi_i, Qj+i, . . . , Qk as constant. Based on section UlI-All we need to solve 
the above problem and prove that the solution is unique. 

Since the numerator and denominator in (fT4l) are concave and affine respectively, (fTSl) is 
a concave fractional program [|7]|. Define a non-negative parameter A, (fT5l) is related to the 
following convex function separating numerator and denominator with help of A. 

F(Q„ X) = k + W\og\l + Gf Q,G,| - A + a,) (16) 

And then define a convex optimization problem as 

r(A) = max F(Qi,A). (17) 

Q,:Q,>0 

We will try to solve (flTI) . and then the solution of (fT5l) can be obtained correspondingly based 
on the following Theorem. 

Theorem 1: The optimum feasible transmit covariance matrix Q* achieves the maximum value 
of dH if and only if F(A*) = F(Q*, A*) = maxF(Qi, A*|Qi > 0) = 0. 

Proof: See Appendix lAl 

Theorem [H gives us insights to solve (fTSl) . We should optimizing (fTTI) at first under a given A 
and then solve the equation ^(A) = to get the optimal A. 
To solve (fTTI) . we can denote 

Gf G, = UD,U^ (18) 

based on the eigenvalue decomposition at first, where Dj G C*^^^ is diagonal with nonnegative 
entries and U G C^^^^^ is unitary. We assume that Dj has L non-zero diagonal entries (1 < 
L < M), which means [Dj]jtfc > for k = 1, . . . , L and [Dj]fcfc = for k = L + 1, . . . , M. 
And then we have the following equation based on I + AB = I -|- BA |l5l: 

log |I -f- Gf Q,G,| = log |I + Q.Gf G,| = log |l + Q.UD,U^| = log |l + U^Q.UD,| 

(19) 

Define = U^QjU. As U is unitary, we have that Tr(Si) = Tr(Qi). Thus, ^ can be 
rewritten as 

G(S„ \) = hi + W\og |I + S,D,| - A (^^^ + a}j (20) 
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As each Sj corresponds to a Qj via the invertible mapping Sj = U^QjU, solving (fTTl ) is 
equivalent to solving the following convex optimization problem. 

Y(X) = max G(Si, \) (21) 

S,;:S,>0 

It is proved in the [5, Appendix II] that the optimal S* to solve (|2TI) is diagonal with [S,*]^^ > 
for A; = 1, . . . , L and [S*]^^ = for A; = L + 1, . . . , M. Thus, ^(8^, A) with diagonal Sj is 

G{S,,X) = k + Wphg{l + [S,]kkWkk)-\ ^SL^ + a}j . (22) 

As (l22l) is concave in Sj, the problem (|2T1) can be solved for a given A by solving the KKT 
optimality conditions, and the solution can be denoted as 

k = l,...,L, (23) 



kk 



_ln(2)A [D,]kk_ 

where = max(x,0). Then the water level A* can be decided by setting Y{X*) = based 
on Theorem \T\ as 



L 



6, + E log 1 + 



ln(2)A* [D, 



' L r ^ 



D, ,fe - A* X !^ ^ + a, = 0. (24) 



A:=l 

Based on (|24l) and (|23l) . the optimal S*^ is derived. Based on the mapping between Sj and Qi, 
finally, the optimal solution of (fT5l) can be derived as 

Q* = USf *U^. (25) 

To prove that the solution of (fT5l) is unique, we only need to prove that A* is unique. We give 
the following Theorem and the proof is given in the Appendix El 

Theorem 2: The derived water level A* in (|24l) is unique and globally optimal for each 
iteration. 

To make the description more clearly, we summarize the energy efficient waterfiUing algorithm 
for optimizing ([H) in TABLE H 

3) Iterative Algorithm: Based on the derivation in section IIII-A2I and the block-coordinate 
ascent algorithm, the energy efficient iterative waterfiUing scheme can be derived as shown in 
TABLE nil and the proof of converge is given as follows. 

Proof of converge: Firstly, during each step, the energy efficient waterfiUing can achieve an 
global maximum EE treating the other users' transmit covariance matrices as constant, the EE 
is non-decreasing with each step. As the EE is bounded, the EE converges to a limit. 
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Secondly, according to Theorem [T] and Theorem |2l the derivation of each step is unique. Based 
on |l33l Sec. 2.7], the set of Qi, . . . , Qk also converge to a limit. □ 

Note that as the proof does not depend on the starting point, we can start the algorithm from 
any starting values of Qi, . . . , Qx- To show the efficiency of the proposed scheme, we give the 
simulation results in Fig. [21 The converge behavior of the proposed scheme is shown and it is 
set that d = 1km, M = A, N = A and K = 10. We can see that our proposed iterative scheme 
converges very fast. It can achieve the optimal EE under nearly five iterations. 



B. Constrained EE Optimization 

After deriving the unconstrained optimal transmit covariances, we can extend it to the general 
constrained case (flOl) with fixed T. Look at the constrained case, the sum capacity is constrained 
and the sum transmit power is limited at the BS. We rewrite them as follows. 

E Tr (Q,) < P^ax 

, K \ (26) 

Cmac ( Hi, . . . , Hj^, E Tr (Qi) 1 > Cmin- 

Define two related problems. The first one is capacity maximization under power constraint. 

(27) 

And the second one is total transmit power minimization under capacity constraint. 



max W log 

{Qi}£i:Q«>0,E£i Tr(Q,,)<Pmax 



i=l 



K 

min E Tr (Q_i) 

i=l 

s.t.t^log I+^EHfQ.H, 

i=l 



(28) 



>C„ 



The first problem (|271) can be solved by the spectral efficient iterative waterfiUing schemes 
efficiently in [5J, while the second power minimizing problem (l28l) can be solved efficiently 



based on |[34l . We denote the solution as Qi and Qi,i = 1,. . . ,K, respectively. To help the 

K _ K _ K ^ ^ 

description, we denote that P* = ^ Tr (Q/), P = E Tr (Q,i) and P^ax = E Tr (Qi). 

2=1 ' 2=1 ' 2=1 V / 

As the numerator of the EE is concave and monotonously increasing and the denominator is 
affine, the EE is quasiconcave, monotonously increasing as a function of the P when P < P* 
and monotonously decreasing as a function of the P when P > P*. We can get the solution of 
the constrained problem as follows. 
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opt 



Q,*,ifP<P* <P^ax 

_ - ,t = l,...,K. (29) 

Q^^iP* < P < Prr 



- 5 

max 



infeasible, ifP > Pmax 

Note that when the solution is infeasible, we will choose Q°''* = Qj,z = 1, . . . as the 
optimal solution in the simulation. This choice has practical significance, as when the rate 
constraint is infeasible, the best choice is to transmit in the available maximum capacity. 

IV. Active Transmit Antenna Selection 

To optimize T, we apply the ATAS procedure in this section to determine it. 

From the standpoint of SE, activating all transmit antennas is always optimal, as higher 
multiplexing and diversity gain can be acquired. However, things change under the energy 
efficient scenario. As multiple antennas consumes higher dynamic power, there exists a tradeoff 
between the power consumption loss and the capacity gain. Thus, ATAS is necessary here. The 
ATAS here is different from the spectral efficient transmit antenna selection, as the conventional 
transmit antenna selection is always utilized in the scenario when the transmit antennas is more 
than the RF chain and the purpose is to employ the selection diversity. Although the ATAS 
can also acquire the selection diversity gain here, its main objective is to compromise the EE 
gain with higher multiplexing/diversity and the EE loss with more transmit independent power 
consumption of multiple antennas. After selecting the active transmit antennas in ATAS, the 
inactive antennas should be turned off through micro-sleep or DTX. 

It is intuitive to see that the exhaust search ATAS is the optimal selection scheme. For each 
possible T C {1, . . . , M}, the exhaust search ATAS calculates the EE based on section Unl and 
then choose the optimal transmit antenna set as follows after comparing the EE. 

Topt = arg max ^(T), 
rc{i,...,Af} 

Finally, the antennas in Topt activated and the inactive ones are switched off . However, the 
complexity of the exhaust search ATAS is too high to implement. Thus, developing low complex 
scheme is important. Let us look at the problem formulation again. As the channel matrices are 
only related to the capacity, the capacity maximizing antenna selection scheme can be borrowed. 
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As shown in ETTl . [|23l . the capacity is highly dependent of the norm of the channel matrices 
and selecting antenna based on norm is highly recommended due to its good performance and 
low complexity. Thus, it is borrowed here, which is shown in Table Hill In the algorithm in 
Table Hill the EE with selected antenna set under different active transmit antenna set size (the 
transmit antenna number) is compared. That is because the power is highly related to the active 
transmit antenna number. For each antenna number, we can choose the channel matrices with 
largest norm. After that, the EE with different antenna number should be compared and then 
determine the best antenna number. 

A. Implementation Issue in Realistic Scenario 

There exists an invisible CSIT problem here. During performing the ATAS, channel matrices of 
all antennas are needed for calculating the Norm and determine the best active antennas. However, 
as the inactive BS antennas should be switched off to save energy, the channel estimation for 
these channels with inactive antennas is impossible. Thus, there might exist time slots in which 
the channel matrices of the inactive BS antennas are not visible at the BS. 

In order to combat this drawback, a possible way is to add one dedicated training period 
to switch on all the BS antennas to obtain the full CSIT of all antennas to help the antenna 
selection. In this case, the power consumption of the training period would decrease the EE. 
Thus, the above definition of EE performs as an upper bound. Nevertheless, there should be 
other low complex schemes to combat the invisible CSIT problem. For example, the statistical 
CSIT can be applied for the ATAS. The performance and cost tradeoff of these schemes should 
be left for the future work. 

V. Simulation Results 

We evaluate the performance under different scenarios to show the effect of parameters. In 
the simulation, pathloss and Rayleigh fading are considered. The parameters are set based on 
H. It is set that W = 5MHz, the noise power is -llOdBm, Pdyn =83W, Psta=45.5W, 77=0.38, 
-Pmax=46dBm and pathloss is 128.1 + 37.6 log^Q d with distance d = 1km (d in kilometers and all 
users are with the same distance). We use "EE w Exh-AS" to denote the optimal energy efficient 
transmission with exhaust ATAS, "EE w Norm-AS" to denote the energy efficient transmission 
with low complexity norm-based successive ATAS, "EE wo AS" to denote the energy efficient 
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transmission with activating all available BS antennas and 'SE' to denote the spectral efficient 
transmission with activating all available BS antennas. Here, schemes with "EE" perform the 
transmit covariances optimization determined according to section |lll] and schemes with "SE" 
perform the transmit covariances optimization based on the capacity maximization |[5| . 

The EE versus pathloss is first evaluated in Figs. [3] and HI where unconstrained case and 
constrained case with Pmax = 46dBm are both considered. We can see that schemes under "EE" 
with ATAS are better than "EE wo AS", and the performance gain comes from the ATAS. In 
Fig. [3l the distance varies from 0.1km to 1km. When the distance is short, e.g. 0.1km, it is 
shown that the schemes with "EE" are all superior to "SE 46dBm", whose performance gain 
comes from the transmit covariance optimization, while the three schemes with maximum power 
constraint perform the same as the corresponding unconstrained ones, whose reason is that the 
globally optimal P* is smaller than 46dBm and thus the constrained and unconstrained scenarios 
can both achieve the global optimality. When the distance gets larger, the "EE wo AS 46dBm" 
degenerates to "SE 46dBm". The reason can be explained as follows. As the distance increases, 
the global optimal energy efficient sum power increases. For "EE wo AS 46dBm", the derived 
P* > 46dBm, and then 46dBm should be applied. 

Look at Fig. |4] then. When the distance is larger than 1km, "EE wo AS" is superior to "EE 
wo AS 46dBm". This is because when the distance becomes larger, the optimal sum power P* 
would be greater than 46dBm, and thus the performance gain comes from the transmit covariance 
optimization. Look at the case when d = 5km, due to the similar reason above, we have that 
"EE w Exh-AS" and "EE w Norm-AS" are better than "EE w Exh-AS 46dBm" and "EE w 
Norm-AS 46dBm", respectively. 

We simulate the effect of capacity constraints in Fig. [5l [6l |7] and the case with M = A, N = 
2, K = 2, d = 1km, Pmax = 46dBm is considered. We can see that the "EE w Exh-AS" always 
achieves the maximum EE, and "EE w Norm-AS" performs very close to "EE w Exh-AS". 
Meanwhile, "EE wo AS" has smaller EE than "EE w Norm-AS", while "SE" has the worst EE. 
The performance gain between "EE" with ATAS and "EE wo AS" comes from the ATAS gain, 
as after selective antennas are determined, turning off the inactive antennas can save the dynamic 
power MaPdyn and then improve the EE. Meanwhile, the gain between "EE wo AS" and "SE" 
shows the efficiency of the energy efficient transmit covariances optimization in section Unl 

The EE gap between different schemes are becoming smaller as the capacity constraint 
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increases. When the capacity constraint is larger than 33 bps/Hz, the four schemes perform 
the same. Look at Fig. [6l the simulated capacity is fixed at 32.2 bps/Hz, and in this case, the 
maximum power is utilized to transmij^. Another observation is that the performances of "EE w 
Norm- AS" and "EE w Exh-AS" both have the multi-stage feature. This feature can be explained 
according to Fig. Ul In order to fulfill the increasing capacity constraint, the active antenna 
number increases, and the active antenna number has a similar multi-stage feature. 

Fig. [8] depicts the EE under different BS antenna configuration, where N = 2, K = 2, Pmax = 
46dBm, Cmin = 0,d = 1km. "EE w Norm-AS" and "EE w Exh-AS" are both monotonously 
increasing as a function of the transmit antenna number M at the BS. The performance gain 
comes from the transmit antenna selection diversity gain with suitable active transmit antenna 
number M^. When the transmit antenna number increases, the probability of choosing channels 
of active antennas with better channel conditions increases. However, look at "SE" and "EE wo 
AS". Their performances degenerate seriously when the antenna number is more than four. The 
reason can be explained by the multiplexing gain of the DPC. There are min(A^ x K, M) spatial 
dimensions for DPC [ l35l and thus the capacity can scales as only N x K for the case with 
M > N X K. For the "SE" with maximum power transmission, the capacity can increase linear 
as the antenna number when M < 4, that is why "SE" has the best EE performance when the 
antenna number is four. Meanwhile, the dynamic power increase linearly with Ma in the power 
part, and then the EE loss with increasing dynamic power is significantly larger than the EE 
gain with capacity increasing when M > 4. Thus, the EE would decrease significantly for "SE" 
and "EE wo AS" then. Furthermore, from the fact that schemes with "EE" and ATAS are much 
superior to "EE wo AS", we can conclude that choosing Ma = M is not optimal for the limited 
user number case here. 

Although it is demonstrated that more antenna number is benefit for the EE with higher 
selection diversity under ATAS, it is worthwhile to note that configuring more antenna would 
cost higher Capital expenditures (CAPEX). In the design of the realistic systems, the tradeoff 
between the EE gain and the CAPEX loss should be taken into account. 

We are interested in discussing the multiuser diversity finally through Fig. [9l which depict 
the EE under different user number, where M = 4, = 2, Pmax = 46dBm, Cmin = 0, d = 1km. 

^In this case, the capacity constraint is infeasible, and the maximum power is employed. 
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We can see that "EE w Norm-AS" and "EE w Exh-AS" degenerate into "EE wo AS", where 
the all four transmit antennas should be active. Moreover, about the multiuser diversity, we can 
see that from Fig. |9]that a similar Mloglog(A^K) scaling law can be acquired. Indeed, in our 
another work ||36l , we analyze the EE scaling law with the help of the Lambert uj function, and 
it is shown that when MaPdyn + -Psta > in this paper, the multiuser diversity of A/^p^'°^|p'^^ 
always holds. We can see that Ma = M is optimal for the large user number case, which is 
distinct from the limited user number case. 

VI. Conclusion 

We propose a novel optimization framework with transmit covariance optimization and ATAS 
to improve the EE in the MIMO BC. Under fixed active transmit antenna set, we transform the EE 
of MIMO BC based on uplink-downlink duality into a concave fractional programming. Based 
on this feature, we propose an energy efficient iterative waterfilling scheme to maximize the EE 
for the MIMO BC according to the block-coordinate ascent algorithm. We prove the converge 
of the proposed scheme and validates it through simulations. After determining the transmit 
covariance under fixed antenna set, we discuss the ATAS to further improve the EE, where 
exhaust search and Norm-based successive schemes are borrowed to select the active antennas 
and switch off the inactive ones. Through simulation results, the effect of system parameters on 
the EE is discussed. 

Appendix A 

As in (fT5l) the numerator is concave and differentiable, and the denominator is convex and 
differentiable. Theorem [U can be directly obtained based on [7, Proposition 6]. 

Appendix B 

As in (fT5l) the numerator is concave and continuous, and the denominator is convex and 



continuous, F{\) is strictly decreasing and continuous based on [|7||. Look at (|24)) . we have that 
F(0) = oo and F(oo) = — oo. Therefore, there exists a unique A* with F(A*) = 0. Furthermore, 
as shown in [7J, in a concave fractional program, any local maximum is a global maximum. 
Therefore, the derived A* is global optimal for each iteration. 
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TABLE I 

Energy Efficient Waterfilling Algorithm 



1) Calculate = I + 77 E Hf Qj Hj , Gi = 

H, I a^I + E Hf Q,H, , a, = ^^^^^ + 

MPdyn + Psta, = H/log|Z,! ; 

2) Define the related parametric convex program in l ll6t 
and {17); 

3) Transform the parametric convex program into diagonal 
forms i20l and l l21t by performing eigenvalue decom- 
position in (|18p : 

4) Solve | |2U by solving the Karush-Kuhn-Tucker optimal- 
ity conditions and obtain the solution S*^ in i23i ; 

5) Calculate the energy efficient water level A* based on 
J24t and determine the optimal S*^ ; 

6) Obtain Q* based on the mapping l |25l > finally; 



TABLE II 

Energy Efficient Iterative waterfilling Scheme 

Initialization: Set Qi = 0, i = 1, . . . , K., 
Repeat: 

Fori = l: K 

1) Calculate Q* based on the energy efficient waterfilling 
algorithm in TABLE HI 

2) Refresh Qi as Q*; 
End 

Until the EE converges. 



22 



TABLE III 
Norm Based ATAS 



Algorithm 

Initialization: Set ^tomp = 0. Sorting the columns of H as 
||H(:,4l))||>...>|iH(:,7r(M))|j. 

For: Ma^l: M 

1) Transmit antenna selection: Choose Ti\f„ — 
{n{l), . . . ,n{Ma)}, and the active channel matrix of 
Ma selected transmit antennas is denoted as Ha/„. 

2) Compute the EE: Calculate the EE based on section 

3) Compare the EE: If ^temp < ^{Ma, Ctcmp = S.{Ma), 
set T = 7ji/„ . 

End For 

Activating the antennas in T and switching off the inactive 




Fig. 1. System model of the MIMO BC. 
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Fig. 2. EE converge behavior of the proposed scheme. 
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Fig. 3. The EE versus distance from 0.1km to 1km, where M = 4,N = 2,K = 2, P^ax = 46dBm, Cmin = 0. 
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Fig. 4. The EE versus distance from 1km to 5km, where M = 4, N = 2, K = 2, P^ax = 46dBm, Cmin = 0. 
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Fig. 5. The EE under different capacity constraints, where M = 4, N = 2, K = 2, Pmax = 46dBm, distance is 1km. 
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Fig. 6. The corresponding capacity under different capacity constraints, where M = A,N = 2,K = 2, Pmax = 46dBm, 
distance is 1km. 
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Fig. 7. Corresponding active transmit antenna number versus capacity constraints. 



February 17, 2012 



DRAFT 



26 




Fig. 8. The EE under different BS antennas, where N = 2,K = 2, Pmax = 46dBm, Cmin = 0, distance is 0.5km. 
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Fig. 9. The EE under different user number, where M = 4, iV = 2, Pmax = 46dBm, Cmin = 0. 
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